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ABSTRACT 

We present results of a study of the stellar age distributions in the sample 
of elliptical-like objects (ELOs) identified at z = in four hydrodynamical, self- 
consistent simulations operating in the context of a concordance cosmological 
model. The simulations show that the formation of most stars in each ELO of 
the sample is a consequence of violent dynamical events, either fast multiclump 
collapse at high zs, or mergers at lower zs. This second way can explain the age 
spread as well as the dynamical peculiarities observed in some ellipticals, but its 
relative weight is never dominant and decreases as the ELO mass at the halo 
scale, M v i r , increases, to such an extent that some recent mergers contributing 
an important fraction to the total ELO mass can possibly contribute only a 
small fraction of new born stars. More massive objects have older means and 
narrower spreads in their stellar age distributions than less massive ones. The 
ELO sample shows also a tight correlation between M v - a and the central stellar 
l.o.s. velocity dispersion, crf^ r - This gives a trend of the means and spreads of 
ELO stellar populations with cr^ that is consistent, even quantitatively, with the 
age effects observationally detected in the stellar populations of elliptical galaxies. 
Therefore, these effects can be explained as the observational manifestation of the 
intrinsic correlations found in the ELO sample between M vir and the properties 
of the stellar age distribution, on the one hand, and M vir and crf^ r , on the other 
hand. These correlations hint, for the first time, at a possible way to reconcile 
age effects in ellipticals, and, particularly, the increase of a/<Fe> ratios with 
a ios"o' w ith the hierarchical clustering paradigm. We briefly discuss the origin of 
the intrinsic correlations shown by ELOs in terms of the adhesion model. 
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1. INTRODUCTION 

Understanding how galaxies form is one of the most challenging open problems in cos- 
mology. Elliptical (E) galaxies are the simplest ones and those that show the most precise 
regularities in their measured properties, with their stellar central l.o.s. velocity dispersions, 
a \osfi ( a mass indicator), strongly correlated with many other of their properties including 
luminosities, sizes, colors and index-line strengths, as recently confirmed by the Sloan Digital 
Sky Survey 'SDSS' (York et al. 2000; Bernardi et al. 2003a,b,c) 

The age distribution of stellar populations in E galaxies is one of the key pieces in the 
puzzle of the origins and evolution of galaxies. It was conventionally thought that these 
populations are coeval systems, formed as the result of a unique burst of star formation 
(SFB) at very high z. Recently, spectral indices have been identified (H^, H 7 , H 5 ) that 
break the age-metallicity degeneracy, allowing for an improved stellar age determination in 
E galaxies through evolutionary synthesis models (see review in Maraston 2003). Even if still 
hampered by uncertainties, these models point now to more massive Es having older mean 
ages (MAs) and lower rates of recent star formation and, also, higher suprasolar a/<Fe> 
ratios than less massive ones (e.g. Jorgensen 1997; Trager et al. 2000b; Thomas, Maraston, 
k Bender 2002; Terlevich & Forbes 2002; Caldwell et al. 2003). The high a/<Fe> ratios 
are actually depressed iron abundances at higher a^ r , rather than a element enhancements 
(Worthey et al. 1992; Trager et al. 2000a), so that the ct/<Fe> ratios are a good measure 
of the timescale for stellar formation (Weiss, Peletier & Matteucci 1995; Thomas, Greggio, 
& Bender 1999; Pagel 2001). The values of these ratios and their correlation with mass 
indicators suggest that an important fraction of the stars in most E galaxies formed on short 
timescales, and that this fraction increases with crf^ r Q . These trends have been confirmed by 
Bernardi et al. (2003c) in their analyses of the SDSS sample of E galaxies, containing to date 
9000 galaxies from different environments. Jimenez et al. (2004) reach similar conclusions 
on the stellar mass - age correlation through a novel statistical analysis of ~ 10 5 galaxies 
from SDSS DR1. 

Two views have historically existed on how Es formed. The modern version of the 
classical monolithic collapse scenario puts the stress on E assembly out of gaseous material 
(that is, with dissipation), either in the form of a unique cloud or of many gaseous clumps, 
but not out of pre-existing stars, with the stellar populations forming at high z and on short 
timescales relative to spirals (Matteucci 2003). The competing hierarchical scenario (e.g. 
Toomre 1977; Kauffmann 1996) propounds that galaxies form hierarchically through succes- 
sive non-dissipative, random mergers of subunits (the so-called galaxy merger tree) over a 
wide redshift range, in such a way that more massive ones (that is, Es) form more likely 
at late time. The age determinations and the interpretation of the a/<Fe> ratios outlined 
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above, as well as their correlations with mass indicators, favor the monolithic collapse sce- 
nario. In fact, the hierarchical scenario tends to predict younger MAs and larger spreads 
in the stellar age distributions of more massive Es (Kauffmann & Chariot 1998; Thomas & 
Kauffmann 1999), see Peebles (2002) and Matteucci (2003) for details and discussions. But 
the monolithic collapse scenario does not recover all the currently available observations on 
Es either. Such are, for example, the wide range in ages their stellar populations span in 
some cases or their kinematical and dynamical peculiarities (Trager et al. 2000a; Menan- 
teau, Abraham & Ellis 2001; de Zeeuw et al. 2002), indicating that an important fraction of 
present-day Es have recently experienced merger events. 

A very convenient approach to reconcile all the current observational background on 
Es within a formation scenario, is to study galaxy assembly from simple physical prin- 
ciples and in connection with the global cosmological model through self-consistent gravo- 
hydrodynamical simulations (Navarro & White 1994; Tissera, Lambas, & Abadi 1997; Thacker 
& Couchman 2000). Individual galaxy-like objects (GLOs) naturally appear as an output of 
these simulations, only star formation processes need further modelling. These simulations 
directly provide the phase space structure and the stellar age distribution of each individual 
GLO at each z, so that the stellar formation rate history of each GLO can be determined. 
This reverses the observational situation where spectral information is available and stellar 
age distributions must be determined through modelling. Also, the structural and dynami- 
cal parameters characterizing each GLO can be estimated. The first step in the program of 
studying the origins of E galaxies through self-consistent simulations, is to make sure that 
they produce elliptical-like object (ELO) samples that have counterparts in the real local 
Universe as far as structure and dynamics is concerned (Saiz, Dommguez-Tenreiro & Serna 
2004, hereafter SDS04). The second step is to show that ELO stellar populations have 
age distributions with the same trends as those inferred from observations and to try to 
understand how these trends arise. These are the issues addressed in this Letter. 



2. SIMULATED ELLIPTICALS 

We have run four hydrodynamical simulations (namely, S14, S16, S17 and S26) in the 
context of a concordance cosmological model, in which the normalization parameter has been 
taken slightly high, <r 8 = 1.18, as compared with the average fluctuations of 2dFGRS or SDSS 
galaxies (Lahav et al. 2002; Tegmark et al. 2003) to mimic an active region of the universe 
(Evrard, Silk & Szalay 1990). Galaxy-like objects of different morphologies form in these 
simulations. ELOs have been identified as those objects having a prominent stellar spheroidal 
component with hardly disks at all. It turns out that the 8 more massive objects identified 
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at z = in S14, S16 and S17, and the 4 more massive in S26, fulfill this condition. We report 
here on age effects in the stellar populations of this ELO sample. This is the same sample 
whose structural and kinematical properties have been analyzed in SDS04, and found to be 
consistent with SDSS data (Bernardi et al. 2003b) for the S16 and S17 subsamples. We refer 
the reader to SDS04 for details on the sample general properties and to Serna, Dommguez- 
Tenreiro & Saiz (2003) for the simulation technique. Star formation (SF) processes have 
been included through a simple parameterization, similar to that first used by Katz (1992), 
that transforms cold locally-collapsing gas at kpc scales, denser than a threshold density, 
Pthrcs, into stars at a rate dp^/clt = c*p gas /£c,, where t g is a characteristic time-scale chosen 
to be equal to the maximum of the local gas-dynamical time, td yn = (47rGp gas )~ 1//2 , and the 
local cooling time; c* is the average star formation efficiency at kpc scales, i.e., the empirical 
Kennicutt-Schmidt law (Kennicutt 1998). To test the effects of SF parameterization, S14, 
S16, and S26 share the same initial conditions and they differ only in the SF parameters (c* 
= 0.1, 0.03 and 0.01, and p th rcs = 6 x 10~ 25 , 1.8 x 10^ 24 ,6 x 10~ 24 g cm' 3 for S14, S16 and 
S26, respectively, that is, SF becomes increasingly more difficult from S14 to S26). To test 
cosmic variance, S17 is identical to S16, except that their initial conditions differ. Supernova 
feedback effects or energy inputs other than gravitational have not been explicitly included 
in these simulations. The role of discrete stellar energy sources at the scales resolved in this 
work is not yet clear. In fact, in the context of the new sequential multi-scale SF scenarios 
(see Elmegreen 2002; Vazquez-Semadeni 2003, and references therein), some authors argue 
that stellar energy releases drive the structure of the ISM only locally at subkiloparsec 
scales. Also, recent MHD simulations of self-regulating SNII heating in the ISM at scales 
< 250 pc (Sarson et al. 2003), indicate that this process produces a Kennicutt-Schmidt-like 
law on average (see also Elmegreen 2002). If this were the case, the Kennicutt-Schmidt law 
implemented in our code would already implicitly account for the effects stellar processes 
have on the scales our code resolves as far as SF rates is concerned, so that our ignorance 
on subkiloparsec scale stellar processes relative to SF (not to metal production or diffusion) 
rates would be contained in the particular values pthrcs and c* take. 

We now briefly comment on the physics of ELO assembly as found in these simulations 
(see Sierra-Glez. de Buitrago et al. 2003). The highly non-linear stages of gravitational 
instability can be described in terms of the adhesion model (Shandarin & Zeldovich 1989; 
Vergassola et al. 1994), based on Burgers' equation (Burgers 1948, 1974), whose solutions 
exhibit time and space scale invariance in a statistical sense in terms of the so-called coales- 
cence length, L c (t), that grows as time elapses, defining the average scale of coalescence in the 
system. Evolution first leads to the formation of a cellular structure that is quasi self-similar 
and not homogeneous, as L c {t) depends on the particular region R within the simulation box 
considered (and, consequently, hereafter it will be written as L c (t, R)). At a given scale, over- 
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dense regions first expand slower than average, then they turn around and collapse through 
fast global compressions, involving the cellular structure elements they enclose and in par- 
ticular nodes connected by filaments, that experience fast head-on fusions. For the massive 
ELOs in the sample this happens between z ~ 6 and z ~ 2.5. These overdense regions act as 
flow convergence regions (FCRs hereafter), whose baryon content defines the particles that 
will end up in a bound configuration forming an ELO at lower z. FCRs contain a hierarchy 
of attraction basins toward which a fraction of the matter flows feeding the clumps they 
host. Another fraction of the matter keeps diffuse. The transformation of gas particles into 
stars at high z mainly occurs through the multiclump collapse ensuing turn around, that 
takes the clumps closer and closer along filaments causing them to merge at very low relative 
angular momentum and, consequently, without orbital delay. This results into strong and 
very fast SFBs at high z that transform most of the available gas in the clumps into stars, 
the exact fraction depending on the values of L c (t, R) at the FCR and of the SF parameters. 
The frequency of head-on mergers decreases after ELO global collapse, notwithstanding that 
at low z mass assembly occurs mainly through merger events. In fact, many ELOs in our 
sample have experienced at least one major merger event at z <1, but a strong SFB occurs 
only if enough gas is still available after the intense thermo-hydrodynamical activity epoch 
at high z. So, our simulations confirm the triggering of SFBs by different dynamical events, 
either at high z (flow singularity formation) or lower z (mainly merger or interaction events 
Tissera 2000). The simulations also show that diffuse gas is heated and expelled at violent 
events involving massive objects, and that a fraction of it is lost to the ELO potential well. 
It forms an X-ray emitting corona whose mass content is similar to the stellar mass content 
bound to massive ELOs (that is, a ~ 50 % of the initial gaseous mass). 

Figure 1 illustrates the star formation rate history (SFRH) of a typical massive ELO in 
the sample (stellar mass M^ Q ar = 3.93 x 10 n M Q at z = 0). Note that at high z this SFRH is 
the sum of many SFBs occurring very close in time. Comparing this SFRH with a detailed 
history of the ELO assembly, it has been found out that in fact the peaks correspond to 
SFBs in different gaseous clumps toward which baryons that form the object at z = flow 
at early times; these clumps merge causing or increasing the intensity of their SFBs, whose 
width is always narrow. Note that the baryon content of individual clumps in an ELO FCR 
is not mixed up before global ELO collapse. Once collapse is over when the universe age 
was t c (t c /t\j ~ 0.3 in this case), further SF occurs through random interaction or merger 
events. An important parameter is then the mass ratio M^ ar ' C0 7M^ ar , giving; the fraction 
of the ELO total stellar mass at z — 0, M^* ar , formed at high z in the ELO collapse event. 
We have found that it increases with the ELO total mass at the halo scale, M vir , or put in 
other words, at fixed p t hres and c*, the fraction of gas at the FCR that is transformed into 
stars at collapse event increases with M vir , leaving a lower gas fraction available to form 
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stars at lower z. At z = the cold gas fraction relative to stars inside ELOs is ~ 2% to 8%, 
depending on their mass. Another important point illustrated by Figure 1 is that most SFBs 
at t > t c are triggered by minor mergers with small gas rich clumps, and not necessarily by 
major merger events, that, nevertheless are very relevant in ELO mass assembly. In fact, the 
ELO in Figure 1 experienced a major merger (mass ratio ~ 0.9) at tf orm /t\j ~ 0.9 or z ~ 0.1 
and, however, no important SFB is apparent in its SFRH, because in this case negligible 
amounts of gas bound to the merging ELOs are available to form stars. We have found that 
major mergers of this kind, that is, involving very gas poor objects, become less frequent as 
z grows, and, in a given z range, as M vir decreases. We have also found that more massive 
ELOs not only form more stars (an issue already analyzed in SDS04), but they also form the 
bulk of their stars earlier on and through stronger SBs than less massive ones. A consequence 
is that the age distribution of stellar populations changes with ELO mass. 

To quantify these differences, the percentiles of ELO stellar age distributions have been 
calculated. For each ELO in the sample, we have calculated the redshifts, Zf, and the universe 
ages, tf, at which the fraction / of the stellar mass at z — 0, M^* ar , was already formed. 
The percentiles at / = 90 can be taken for a measure of the amount of low z star formation; 
the difference At = t 75 — t w as an estimation of the global width or timescale for ELO star 
formation; t 50 is an estimation of the MA of the population. For any / a trend exists with 
M vir . The observational age effects with crf^ arise because M vir and af^ T are on their turn 
tightly correlated (we have found for the ELO sample ogj, = (0.057 ± 0.029)M v °; r 28±0 - 02 ), so 
that cfosjo * s an empirical virial mass estimator, see Saiz (2003); Saiz et al. (in preparation). 
Moreover, a^o an d the ELO stellar mass content M^* ar are also closely correlated (see details 
in SDS04), so that a trend of tf and Zf with M^* ar also exists. As an illustration of these 
trends, in Figure 2 we plot t^o and At versus o~f^ r for the ELO sample. We see that more 
massive ELOs have older MAs and narrower spreads in the distributions of their stellar 
populations. Figure 2a compares adequately well with relative MA determinations through 
population synthesis modelling for Es in different <J^ r bins. For example, for the isophotal 
populations in the E sample from Caldwell et al. (2003, their Tables 1 and 9), the MA 
averages for Es with af^ in the ranges (260 - 158) and (158 - 100) km s -1 differ in ~ O.ltu, 
consistent with that found for the global populations of the ELO sample 1 . Figure 2b shows 
that width determinations from a/<Fe> ratios are consistent with ELO widths except for 
the most massive ones, at, say, of*^ ~ 240 km s _1 (crosses in Figure 2b are from Thomas, 
Maraston, & Bender (2002, Figure 6), maybe because the most massive ELOs have been 
dynamically shaped at low z through major merger events involving smaller ELOs of ~ half 



1 Trager et al. (2000b) get a difference of ~ 0.2 — 0.3tu for the central stellar populations of the Gonzalez 
(1993) E sample. 
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their final mass, so that the stellar populations of these massive ELOs at z = reflect the 
properties of the stellar populations corresponding to the smaller merging ELOs. Note that 
the trends in Figures 2 appear in subsamples with very different SF parameterizations, so 
that the trends are independent of the particular details of the SF implementation. 

3. DISCUSSION 

The simulations we report on in this Letter not only show that ELO stellar populations 
have age distributions with the same trends as those inferred from observations, but they 
also provide us with clues on how these trends arise. They indicate that ELOs are assembled 
out of the mass elements that at high z are enclosed by those overdense regions R whose 
local coalescence length L c (t, R) grows much faster than average and whose mass scale (total 
mass enclosed by R, Mr) is of the order of an E galaxy total (i.e., including its halo) mass. 
The virial mass of the ELO at low z, M vir , is the sum of the masses of the particles that 
belong to R and are involved into the ELO merger tree. Note that this sum is determined 
by the local realization of the initial power spectrum of density fluctuations at R, that, on 
its turn, also determines the total energy of the particles in the ELO, E (observationally 
seen as velocity dispersion af*^), and the growth rate of L c (t,R). As a consequence, when 
L c (t,R) grows faster, M vir and E are higher (see also Saiz et al., in preparation, and refer- 
ences quoted therein). Now, when L c (t,R) grows faster than average at a given region R, 
the local time unit at R is shorter than at other R's where it grows slower. And so, the 
thermo-hydrodynamical activity at high z at R is more intense than at other R's and, in 
particular more stars form and gas becomes exhausted earlier on, making it difficult further 
star formation at low z merger events. This is how the intrinsic correlations age distribu- 
tion - M vir - E arise in an ELO sample, whose observational counterpart are likely the age 
distribution - a^ correlations shown by elliptical samples, as discussed in §1 and §2. 

The simulations we analyze also suggest that a fraction of the stars in E galaxies could 
have been formed at lower redshifts. A related interesting point is the possibility that mass 
assembly and dynamical shaping at low z is not necessarily accompanied by strong SFB ac- 
tivity. For example, the two more massive ELOs in S14 and S16, or the most massive in S17, 
have been structurally and dynamically shaped at a low z (z ~ 0.1 and 0.65, respectively), 
however only very modest starbursts resulted from these dynamically very violent events, 
because of gas consumption at high z. 

Some interpretations of the hierarchical clustering paradigm consider that the [a/ <Fe>] 
- crj^gQ correlation is not consistent with this paradigm. The results reported in this Letter 
indicate that they are indeed, provided that the bulk of the stellar populations forming ELOs 
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have been formed in different subunits at high z and merged together very soon after (see 
Thomas, Greggio, & Bender 1999; Sierra-Glez. de Buitrago et al. 2003). In this Letter we 
show that this fast clumpy collapse follows from simple physical principles in the context of 
the current ACDM scenario. 
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ported by the MCyT (Spain) through grants AYA-0973, AYA-07468-C03-02 and AYA-07468- 
C03-03 from the PNAyA. We thank the Centro de Computacion Cientffica (UAM, Spain) 
for computing facilities. AS thanks FEDER financial support from UE. 
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Fig. 1. — Star formation rate history of the more massive ELO formed in S16 as a function 
of universe age in units of the current universe age 
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Fig. 2. — (a) Age of the universe in units of the actual universe age at which the 50 per 
cent of the total ELO stellar mass at z = was already formed, versus their corresponding 
stellar central l.o.s. velocity dispersion. Filled triangles and pentagons stand for S16 and 
S17 ELOs; open squares and circles for S14 and S26 ELOs, respectively, (b) Same as (a) 
for the width of the stellar population age distribution. Crosses are width estimations from 
elliptical data, see text. 



